Exact diagonalisation study of charge order in the quarter-filled 
two-leg ladder system NaV205 
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The charge ordering transition in the layer compound NaV205 is studied by means of exact diago- 
nalization methods for finite systems. The 2-leg ladders of the V- Trellis lattice are associated with 
one spin variable of the vanadium 3d-electron in the rung and a pseudospin variable that describes 
its positional degree of freedom. The charge ordering (CO) due to intersite Coulomb interactions 
is described by an effective Ising-like Hamiltonian for the pseudo-spins that are coupled to the spin 
fluctuations along the ladder. We employ a Lanczos algortihm on 2D lattice to compute charge 
(pseudo-spin) and spin-correlation functions and the energies of the low lying excited states. A 
CO-phase diagram is constructed and the effect of intra-ladder exchange on the CO transition is 
studied. It is shown that a phase with no-longe range order (no-LRO) exists between the in-line and 
zig-zag ordered structures. We provide a finite-size scaling analysis for the spin excitation gap and 
also discuss the type of excitations. In addition we studied the effect of bond-alternation of spin 
exchange and derived a scaling form for the spin gap in terms of the dimerization parameter. 
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I. INTRODUCTION 

The interplay of charge ordering (CO) and spin struc- 
ture and dynamics has recently been in the focus of in- 
terest especially in the layered 3d-perovskites like cubic 
and bilayer manganites. A CO transition is caused by 
inter-site Coulomb interactions dominating the kinetic 
hopping energy of the electrons of transition elements be- 
low a critical temperature Tc. It localizes the electrons 
and thus defines an associated low temperature spin lat- 
tice which then may exhibit magnetic ordering and pos- 
sibly low dimensional spin excitations. Such an interplay 
of charge order and spin excitations, supplemented by 
a lattice distortion is present in the Trellis lattice com- 
pound NaV2 05 . Its basic building blocks consist of two 
leg ladders whose rungs are occupied by one d-electron 
on the average (Fig.^). It is an insulating compound 
with an effective charge transfer gap given by the hop- 
ping ta across the rung. The electronic structure of this 
compound was investigated in Refs.( |]l|- ||]) in detail 
and hopping parameters for a tight binding model were 
derived. Susceptibility measurements Q NMR- experi- 
ments j^] and x-ray analysis have shown that a spin 
gap evolves below Tc ~ 33°-?^ at the same time when 
NaV2 05 undergoes a CO- transition and lattice distor- 
tion. Sofar there is no complete understanding of both 
the low temperature structure and the origin of the spin 
gap. 

A model for the coupling of CO to spin degrees of free- 
dom in NaV2 05 was introduced in Rcf.( Q) based on 
a pseudospin representation for the rungs such that the 



CO transition can be described within a generalized 2D 
Ising model. Many proposals for the proper low temper- 
ature structure have been made. According to the model 
by Liidecke et al |^ only every second ladder has zig-zag 
CO whereas the other ladder stays disordered. This CO 
is also supported by neutron scattering results jl^ dis- 
cussed in Ref.( |l|]). On the other hand, this charge dis- 
tribution disagrees with an x-ray anomalous scattering 
measurement which indicates charge modulation along 
the b-axis [^. Furthermore the structure in Ref.( 
leads to three inequivalent V- sites in valence states , 
V^'^^ and which is incompatible with NMR results 
[Q which show that only two inequivalent V-sites exist. 
A solution of this puzzle has been proposed in Ref.( [|3|) 
where it was shown that this structure actually leads 
to asymmetric and incomplete charge order where two 
of the inequivalent sites have almost identical valences. 
However incomplete charge order complicates the naive 
application of conventional exchange models for the spin 
dynamics. Therefore it is first necessary to study the 
influence of charge ordering on the superexchange mech- 
anism. It was shown [ p^ how charge correlations in the 
rungs and CO progressively reduce the superexchange of 
spins which occupy rung molecular orbitals. For this cal- 
culation to be carried through, some mean field like ap- 
proximations in the charge (pseudo- spin) variables have 
to be made to start from a restricted orbital space for 
the calculation of exchange integrals. Finally we mention 
that some recent results jl^ emphasize the fully zig-zag 
CO in the ladders below Tc for NaV205 , but the stack- 
ing perpendicular to the ladders for this model is still 
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controversial. 

As long as the critical temperature (Tc) is finite one 
may guess that only classical fluctuations are responsible 
for the phase transition phenomena. However an esti- 
mate of the critical point in the whole phase diagram will 
show whether the quantum fluctuations should be con- 
sidered. Since the scale of hopping energies and Coulomb 
interactions are in the order of 1 eV then the critical tem- 
perature in this unit is Tc=0.0028 eV. This shows that 
our model is very close to quantum critical point and 
quantum fluctuations are very important in the phase 
transition. Moreover the results at T=0 can explain the 
phase transition phenomena to a very good extent. It 
is therefore desirable to do a fully quantum mechanical 
calculation which treats the charge (pseudo- spin t) and 
spin (T variables of the V-V rungs on the same level. Such 
an approach is presented in this work and carried through 
by using exact diagonalisation methods with the Lanczos 
algorithm. In Sec. II we introduce the pseudospin model 
that describes the charge ordering transition and the cou- 
pling of intra-rung charge fluctuations to the spins. Sec. 
Ill is devoted to the description of the numerical methods 
employed. Sec. IV gives the computational results for the 
correlation functions and the conjectured phase diagram 
of CO, furthermore the influence of the spin exchange 
on the critical CO parameter is studied. In addition we 
investigate the finite size scaling of the spin excitation 
gap and its dependence on the degree of charge order. 
Moreover we have studied the effect of spin exchange 
bond-alternation on the energy gap and the scaling of 
spin gap for small value of dimerization 7. A summary 
of our results is given in Sec. V. 



II. CHARGE ORDERING AND PSEUDO-SPIN 
MODELS 

An attractive feature of the NaV205 compound that 
simplifies its theoretical study, is the fact that the charge 
ordering happens to be in an already insulating mixed 
valent (MV) state. It is characterized by a nonintcger 
valence 4.5-1- of each vanadium atom of the rung which 
therefore is occupied by only one d- electron on the av- 
erage. The ordering is due to the fact that the large 
intra-rung hopping t^ leads to a half- filled (rung-) bond- 
ing band that is separated into an (empty) upper and 
(filled) lower Hubbard band. The effective insulating gap 
is then a charge transfer gap approximately given by t^. 
Although this scenario is suggestive, a detailed theory for 
the insulating state of NaV205 is still lacking. Due to the 
charge gap, each rung accommodates only one d-electron 
which can occupy the right (R) or left (L) position of the 
rung. As first proposed in Ref . ( |^ ) , this charge degree of 
freedom can be described by an Ising pseudospin variable 
Tz and its resonating behaviour by the transverse isospin 
variable Tx- The latter is responsible for the tendency 



towards the formation of a bonding state. The Coulomb 
interactions between d-electrons on different rungs may 
then be mapped onto Ising-like interactions on the one- 
particle subspace of the rungs. Altogether the charge 
degrees of freedom are then described by 
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FIG. 1. Schematic representation of the crystal structure 
of a vanadium layer of NaV205 . The vanadium ions are de- 
noted by filled circles. The couplings in Eqs.(^ & ^ are also 
shown: the hopping amplitude (ta) on each rung, effective 
pesudo-spin K and spin J couplings on each ladder and be- 
tween ladders K',J'. Different configurations that we have 
considered in exact diagonalization are denoted by: (a) A 
single ladder is considered with A'' — 28. (b) Two interacting 
ladders via K',J' compose a lattice of N — 24. (c) Three 
interacting ladders and in (d) four interacting ladders. 

The first term describes the intra-rung hopping, which 
is the largest hopping element to, see Ref.( [Q). The 
terms ~ K, K' are effective intra- and inter-ladder 
Coulomb interactions, respectively |^ (Fig.|l|), which con- 
tain also the effect of virtual inter-rung hopping. An in- 
teresting particular case of this Hamiltonian occurs if we 
consider only a single ladder (i.e. set K' = 0). Then this 
Hamiltonian describes the ID Ising model in a transverse 
field of strength t^ whose exact solution, including the 
correlation functions (T^r|) and (t^t|) (j=i±l), is known 
|16|. It has been discussed in the present context in Ref.( 
|14|). The model is able to describe the appearance of 
the zig-zag charge ordering along the ladder as a func- 
tion of the control parameter K/ta; strictly speaking an 
infinitesimal staggered field due to inter-ladder coupling 
is necessary to stabilize the order parameter |Q. Due 
to the peculiar structure of the Trellis lattice (rungs of 
neighbouring ladders are shifted by half a lattice constant 
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b) the inter-ladder effective Coulomb interaction K' leads 
to a frustration in the zig-zag CO state. If it exceeds a 
critical value K^, then the zig-zag CO is destroyed and in- 
line CO would be prefered. This transition corresponds 
to a change in the modulation vector of the pseudo-spin 
density t^{Q) from Q = , to Q = 0. The Hamilto- 
nian of Eq.(|^) with K>0 describes a 2D frustrated Ising 
model in a transverse field (ITF) on a trigonal lattice. 

In this work our aim is to present a more detailed study 
of the correlation functions and phase diagram of this 
model. 

It was shown in Ref.( ||l^) that the spin exchange be- 
tween the ladders is strongly influenced by the charge 
fluctuations and the charge order parameter. In this in- 
vestigation the intra-rung charge fluctuations have been 
eliminated by using the ID- Ising model exact solution. 
This is only possible for one ladder and not in the 2D 
frustrated Ising model appropriate for the Trellis lattice. 
Therefore we will keep the general interaction term of 
charge and spin degrees of freedom as proposed in Ref . ( 
@) 

Hcs^J ^[(l + TD(l + ri)aV^' 

+ (l-r:)(l-Ti)aV^] 

+ (l-T^)(l+Ti)aV^] (2) 

The first and second term describe the coupling of 
charge fluctuations to the intra- and inter-ladder ex- 
change, respectively. For clarity we use explicitely the 
R,L pseudospin projectors i(l ± Tz) in the Hamilto- 
nian. Approximate expressions for the superexchange 
constants J, J' have been given in Ref. ( ^ ) . In gen- 
eral we ignore superexchange paths along the ladder di- 
agonals, although this is not completely justified ji); we 
will consider its effects in a few cases. So far we have 
introduced a large number of coupling constants as free 
parameters in these models. It is not our purpose to com- 
pute the properties of them for any range of the param- 
eters involved. Instead, we shall be interested in those 
values that are physically relevant for the underlying 
NaV205 compound. They were obtained in LDA+U cal- 
culations Q assuming full charge order. In this method, 
the 3d-occupation of each vanadium site is fixed and an 
empirical on- site Coulomb interaction U is added to the 
LDA potential which shifts occupied 3d-states downward 
and unoccupied 3d-states upward with respect to the 
Fermi level and in this way leads to a proper insulat- 
ing state for NaV205 . From these LDA-I-U calculations 
Q the following hopping and exchange constants were 
obtained: ta= -340 meV, J = 30 meV and J' = -5 meV, 
i.e. there is antiferromagnetic (AF) intra-ladder and fer- 



romagnetic (FM) inter-ladder exchange. As the effective 
inter-site Coulomb constants are not known, they are left 
as variable parameters in a range of up to 1 eV. The on- 
site Coulomb interaction, which does no longer appear in 
our effective low energy Hamiltonian, has the magnitude 
t/ = 4 eV. 

III. LADDER HILBERT SPACE AND 
COMPUTATIONAL TECHNIQUE 

In our numerical computation, each rung i is associ- 
ated with the four dimensional Hilbert space {Ici) €3 
where ai,Ti =t, |. We classify the states used according 
to their total spin quantum number az and restrict to 
the space with (Tz—O for the ground state calculations. 
Note that the total pseudo spin z-component t^— is 
not conserved due to the first term in Eq. (|l|) . As a conse- 
quence of this, all states in the pseudo-spin sector have to 
be included in the calculation. Therefore the maximum 
number of sites which we treat is N=28. They may be 
arranged in different structures, some of them are shown 
in Fig.|l|. In Fig.(|l|-a) we have plotted an isolated single 
ladder where N=28. The intra-ladder coupling constant 
K (which is the result of Coulomb interaction on neigh- 
bouring rungs) exists on each leg of the ladder. Similarly 
the exchange constant J on each leg represents the cou- 
pling of the spin degrees of freedom of each electron. The 
hopping parameter ta is defined on each rung. Most of 
our calculation on the single ladder is considered with 
periodic boundary condition (BC) in the legs direction. 
We have only considered the open BC when we compare 
the correlation functions of a single ladder and two cou- 
pled ladders. The interactions of two ladders is assigned 
in Fig.(^b) with inter-ladder couplings K' and J' which 
is denoted by dashed lines. All of results presented here 
on two coupled ladders are in open BC. The interaction 
between three and four ladders is depicted in Fig.(|^c 
and d). The results of computations on this latter con- 
figurations are not presented here because the lattice di- 
mension on each direction is very small and we cannot 
avoid boundary effects. However this results help us to 
see what will be the effect of more ladder interation and 
find a qualtitative phase diagram for our model. 
We use a single ladder to study the finite size scaling 
of various quantities as function of the number of sites 
N. We also use some more two dimensional structures as 
shown in Fig jl] corresponding to various N to investigate 
the effects of inter-rung interactions, especially concern- 
ing the phase diagram as a function of K,K'. We rep- 
resent correlations for open boundary conditions on the 
ladder since for the former the correlation functions are 
symmetric around the midpoint of the ladder and there- 
fore presumably contain less information. The periodic 
boundary condition is applied in the b-direction and used 
to compute the energy levels. 
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IV. RESULTS OF EXACT DIAGONALIZATION 

In this section we present the results of our calcula- 
tions for correlation functions of the Hamiltonian (H) 
with both spin and pseudo-spin degrees of freedom [H — 
He + Has)- From this information, we are able to con- 
jecture a phase diagram in the {K, K')- plane. Further- 
more, we also study the influence of the exchange on the 
critical value Kc for charge ordering in the single ladder. 
Subsequently, we come to the important question about 
the nature of excited states (charge type vs. spin type) 
and discuss the finite size scaling of the energy gap in a 
single ladder. We have also considered the effect of bond- 
alternation (dimerization, 5) of the spin exchange on the 
single-ladder system. This shows the scaling behaviour 
of spin gap due to 7. 
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FIG. 2. Pseudospin correlation function ({r^(0)r^(r))) for 
a sine le ladder (Fig.(|l]-a)) versus r, the distance between 
rungs. Periodic boundary condition are considered along legs 
of ladder so that the correlations are symmetric around the 
mid-point. By increasing K, an evolution from a disordered 
phase (no-long-range order) at small value of K (< 0.35) to 
long-range order happens for large value of K (J=30 me V) . 

A. Correlation functions and CO transition 

The pseudo-spin interaction K > along the lad- 
der legs is of "antiferro" type because physically the d- 
electrons sitting on the same side of adjacent rungs have 
a larger repulsion than if they were located across the 
ladder diagonal. Thus the second term in Eq.(|l|) tends 
to localize the electrons on opposite sites of the ladder in 
zig-zag fashion (Fig.^. On the other hand the first term 
in Eq. (|l|) which describes the hopping energy within the 
rung prefers the disordered mixed valent state where the 
d-electron occupies the symmetric rung-bonding state 
with equal probability of sitting on the L and R sites. 



This state is realized for K = (and ta as given above 
from LDA-I-U calculations) . If K increases above a crit- 
ical value Kc — ta, eventually the ordered state with a 
zig-zag correlation is preferred. This behaviour is clearly 
exhibited in Fig.^ for the single ladder (K'=0). 

The correlation function for a single ladder with peri- 
odic BC is plotted for N=28 in Fig.(^). For small values 
of K (K< 0.4) there is no LRO and correlations decay 
very fast. By increasing K (K> 0.4), the correlation 
function shows zig-zag (alternating) LRO where the am- 
plitude of the correlations is approximately constant for 
all distances (r). The critical interaction can be seen 
to be Kc= 0.35-0.4 eV ~ ta. The coupling to adjacent 
ladders (third term in Eq.(|^)) would however prefer the 
"ferro" type arrangement, which leads to a frustration of 
charge (pseudo-spin) order. The effect of this frustration 
is clearly visible in Fig.^ and Fig.^ for the correlation 
functions with K' >0, most directly for the 2"'' neigh- 
bour rung (r=2). Its correlation for K' = is of ferro 
type, then increasing K' becomes first antiferro-type and 
finally, when K' is large enough so that the in-line or- 
dering is prefered, it becomes of ferro type again. This 
non-monotonic behaviour of (r^(i)T^(i-|-2)) is due to the 
competition of different interaction paths, either directly 
along the legs of the given ladder, involving K, or indi- 
rectly via adjacent ladders which involves K'. 

A similar behaviour is seen in the two-dimensional lat- 
tice by considering three or four ladders interacting via 
K'. In this case the number of rungs is reduced to four 
and three respectively in each ladder (Fig.|l]-(c) and (d)). 
We shall explain the case of the two-ladder system in 
more detail. In this case the system is composed of two 
6-rung ladders with open boundary conditions imposed. 
The pseudo-spin correlation functions are plotted in Fig.^ 
and Fig.^. We have plotted (r^(O)r^(r)) along the b- 
direction for both single and double ladder systems. The 
case of a single-ladder is added for comparison. In Fig||, 
the value of K is fixed to K=0.35 < Kc. The case of 
single-ladder (K'=0) shows a rapid decay with distance 
which verifies that no LRO exists in the ladder. By turn- 
ing on the K' coupling, the second ladder interacts with 
the first one. We have plotted the correlation functions at 
the same points of the single-ladder for different value of 
K'. When K'=0.175 or 0.35 eV, the qualitative behaviour 
does not change which verifies the absence of LRO. By 
increasing K' the ordering of the first neighbour (r=l) 
goes from anti-ferro to ferro type. For the 2"'' neighbour 
(r=:2) this ordering goes from ferro- to anti-ferro and fi- 
nally again to ferro- type for enough high value of K' (in 
this case K' > 1.25). As mentioned above, the compe- 
tition between different paths of interactions via K and 
K' is responsible for this behaviour. It happens for 0.35 
< K' < 0.7 eV which is a phase with no LRO and will 
appear in the phase diagram (Fig.^). For high value of 
K' an in-line LRO appears on each ladder of the model. 
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FIG. 3. Pseudospin correlation function (('r^(0)r^(r))) for 
both single (Fig.([l|-a)) and double ladder (Fig.(|ll-b)) versus 
distance (r). In both cases the intra-ladder coupling (K) is 
fixed to K=0.35 < Kc- Then the correlations decay very fast 
and do not show long-range zig-zag ordering. By increasing 
K' (inter-ladder coupling) deviation from single ladder results 
start and leads to a disorderd regime for an intermediate value 
of K'=0.35. For sufficient size of K' an in-line ordering will 
appear. (Correlations computed with open boundary condi- 
tion, J=30 meV and J'=-5 meV). 

In addition, Fig.^ shows a similar behaviour like FigJ|, 
while K is fixed to K=0.5 eV >Kc in the zig-zag ordered 
phase of a single-ladder system. The correlation function 
shows that for the single-ladder (K'=0), a true zig-zag 
LRO exists. As before, a small value of K' does not 
change the ordering of the single-ladder (K'=0.25 and 
0.5 eV). For an intermediate value of K', the disordered 
(no LRO) phase wiU appear (0.5 < K' < 1.5 eV). Then 
by increasing K' a complete in-line ordering starts to sta- 
bilise. 



B. Phase Diagram for CO 

As we have seen in the previous sections, the tran- 
sition from antiferromagnetic (zig-zag) to ferromagnetic 
(in-line) type CO depends on the ratio of K to ta and 
on the "frustration" ratio K'/K. As the 2"^ neighbour 
(r=2) correlation along the ladder shows the transition 
is not abrupt but gradual. This signals an intermediate 
region in the K, K' plane which has to be considered as a 
phase of no LRO. We believe that this phase is very sim- 
ilar to the "floating phase" discovered in the 2D-ANNNI 
(anisotropic next nearest-neighbour Ising) model on a 
square lattice even though in this case it exists only at fi- 
nite temperature (e.g. Ref.( 0)). The inherent geomet- 
ric frustration of the Trellis lattice might stabilize it even 
for zero temperature. The phase diagram obtained from 



systematic study of the correlation function is shown in 
Fig.^. For simplicity, let us neglect the spin coupling and 
consider Eq.(|l|) in two marginal cases. First, if K' — 0, 
Eq.(|l]) reduces to a one-dimensional ITF model. Then 
for K<Kc there is no LRO for charge degrees of freedom 
while a 2"''^-order quantum phase transition will occur at 
Kc to an anti- ferroelectric (zig-zag) LRO. In the other 
extreme case where K=0, Eq.(|l|) is reduced to a two- 
dimensional ITF model on a square lattice. Obviously for 
small value of K', the first term in Eq.(^ prevents the ex- 
istence of LRO for charge degrees of freedom. If K' goes 
to large values respect to ta, the two-dimensional square 
lattice shows a ferro-electric (in-line) ordering. Therefore 
the phase diagram of Eq. ([|) in the K-K' plane consists of 
two boundary lines which divide a phase with no LRO in 
the intermediate regime from the other in-line and zig-zag 
LRO . This phase boundary is an extension of two critical 
points Kc and K'c in the K-K' plane which is supported by 
numerical data on two-dimensional clusters (Figs.^ and 
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FIG. 4. This plot is similar to Fig.(^ but for different val- 
ues of K and K'. Here K is fixed to K=0.5 > Kc which shows 
long-range zig-zag order on the single ladder (K'=0). By turn- 
ing on the inter-ladder coupling (K'), no qualitative changes 
are observed for small value of K'=0.25. By increasing K', 
the system goes through a disordered phase for intermediate 
values of K'=1.0 and then enters to the in-line ordered phase 
for large value of K'=1.5. (Correlations are computed with 
open boundary condition J=30 meV and J' =-5 meV). 

The phase boundary for the upper in-line phase and 
the lower zig-zag phase are obtained by requiring a 
threshold value of 90% for the correlation of each type. 
As critical values we obtain Kc ~ 0.39 eV (the effect of 
spin coupling has been considered) and K^ is the critical 
value of a two-dimensional ITF model on a square lattice 
including the spin coupling by Eq.(|^). If we neglect the 
effect of Eq.(|) then K^ -0.054 eV These phase 
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boundaries enclose the disordered phase with no clearly 
discernible long range correlation. Although specifying 
the properties of the intermediate no LRO phase is be- 
yond our numerical computations our data show that the 
period of oscillation of correlation functions has doubled 
in comparison to the zig-zag ordered phase. This can be 
seen in Fig.(|) at K'=0.7 and in Fig.(|) at K'=1.0. 



K' 



in-line ordering 



K' 



^ zig-zag ordering 



Kc K 

FIG. 5. Charge-order (CO) phase diagram in the space 
of intra-(K) and inter-ladder (K') couplings. Kc repre- 
sent the critical value of K in a single ladder configuration 
(Fig.(^a)). On the vertical axis Kj. is the critical value of a 
two-dimensional (ITF -I- spin exchange effect) on a square lat- 
tice. The dashsed lines show the qualitative phase boundary 
between three different phases where zig-zag ordering happens 
in the lower-right part and in-line ordering in the upper-left. 
A phase with no LRO separates the two ordered phases. 

The phase boundaries for CO also depend on the size 
of the exchange constants. This has been studied for 
the point (K,K') = (Kc,0), corresponding to the single lad- 
der case in more detail as shown in Fig.^ The critical 
value Kc for development of long range order is defined 
by the simultaneous crossing point where curves for dif- 
ferent ladder lengths (N=20-28) meet. One can clearly 
see that Kc shifts as a function of the exchange J along 
the ladder leg. This can be explained as follows. In the 
disordered phase there is a certain amount of stabilisa- 
tion energy of the spin singlet. The incipient zig-zag CO 
which localizes the d-electrons on opposite sites of the 
ladder diagonal tends to reduce this stabilisation energy 
and hence, if J increases the critical Kc for CO to set in 
also has to increase. It is shown for two different values 
of J (30, 100 meV) in Fig.|. In the case of single-ladder 
by setting J = 0, we have a one-dimensional ITF model 
where Kc=ta. It is shown in Fig.^a that the intersec- 
tion of the pseudo-spin correlation function for different 
values of N (20, 24, 28) is at Kc=0.39 eV where J=0.03 
eV (ta=:0.34 eV). A similar curve in Fig.|^b shows that 



the intersection is at Kc=0.40 eV where J is increased to 
J=0.1 eV at fixed to. This verifies the monotonic depen- 
dence of Kc on the exchange constant J . 

< t'(1) t '(R) > al R=N/2 + 1 
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FIG. 6. The pseudo-spin correlation amplitude at fixed dis- 
tance R=-y versus the intra-ladder coupling (K) for single 
ladder. Different sizes of ladder (N=20, 24, 28) are imple- 
mented to track the crossing point of curves, which is the 
transition point, (a) The effect of spin exchange (J) is consid- 
ered with the value of J=30 meV then the crossing point hap- 
pens at Kc ~0.39. (b) In this case J=100 meV which causes 
an increase to the value of Kc ~=0.4. (Periodic boundary 
conditions are considered along the legs of ladder) . 



C. Nature of excited states, finite size scaling and 
effective exchange 

We have implemented the Lanczos method in differ- 
ent Sl^i subspaces to determine the structure of excited 
states. For the single-ladder system, the ground state 
and first excited state are always spin singlet and triplet 
state respectively, as long as we do not include the ex- 
change path along the ladder diagonals. This is true even 
in the coupled ladder system (2D Trellis lattice). As we 
have a Heisenberg exchange term, the singlet- triplet gap 
has to vanish in the thermodynamic limit N— > cxd. 

However, from the finite size scaling behaviour we may 
extract the "effective exchange" along the ladder for a 
pure spin model, which is modified from its bare value 
J due to the effect of intra-rung charge fluctuations. In 
Fig.^a the scaling behaviour of the singlet- triplet gap 
A for a single ladder is shown as a function of 1/N 
(N=16, 20, 24, 28) for various interaction paramters K. 
All curves show a decrease in the energy gap with increas- 
ing N, which leads to gapless behaviour by extrapolation 
to oo. We use a scaling form for the energy gap 
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jiven by 



(3) 



where Aqo ^0 for 'N—>- oo, a is constant and v a scaling 
exponent. These values are given in Table 1 for various 
K and we see that < 1 for K < Kc and > 1 for K > 



K 


a 




0.3 


0.03379 


0.75224 


0.35 


0.03332 


0.88581 


0.392 


0.03633 


1.04005 


0.4 


0.03683 


1.06949 


0.45 


0.03699 


1.21717 


0.5 


0.02862 


1.23286 



Table 1. Numerical values of the coefficients in the scal- 
ing function of the energy gap in Fig.^a. The scaling 
exponent is < 1 for K < Kc — 0.39 and > 1 for 
K > Kc- The scaling form is A = Aqo + (Data has 
been obtained by imposing a least square root error of 
10~^ for the scaling function ). 
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FIG. 7. The energy gap to the excited state versus 
1/N of the single ladder (Fig.(^a)) for different values of in- 
tra-ladder coupling (K). (a) All data show gapless behaviour 
in the A'' — > oo. The spin exchange is J=30 meV . (b) In ad- 
dition to J we have added diagonal exchange terms (Jd= 32.7 
meV) on each plaquette of ladder. For K<Kc=0.39, we ob- 
serve qualitative behaviour similar to part (a) where the first 
excitation is a spin triplet. But for K>Kc the behaviour is 
completely different because the excitations are charge type. 
The excited state has the same total spin Stot=0. 

This shows the different curvature of A versus 1/A^. 
(To have the best fit of the data in Figj^a at fixed ly — 1 



we have considered more terms in the curve fitting form; 
A{N) = Aoo + ^ + +_% which is justified for the 
case of pure spin- i chain [^). For a pure spin-i system 
without coupling to charge degrees of freedom, one has 
iy~l and the coefficient Ci is related to 



Ci = 2n^ Jeff 



(4) 



Comparison with our Ci coefficient from the scaling plots 
in Figj^a gives us the variation of the effective exchange 
constant Jeff(K) with the interaction K that controls the 
charge order (for a given 1^). This variation is shown 
in Fig.|^. A steep drop in Jeyj(K) is seen in the region 
where the charge order sets in (K~0.39). Eventually for 
K^ta when zig-zag CO is almost complete, Jeff has to 
approach zero because in the zig-zag structure the ex- 
change J along the legs is ineffective. 

The situation changes when we also include exchange 
along the ladder diagonals. In this case we will add 
an exchange term similar to the first term of Eq.(H) 
along the ladder diagonal in each plaquette (not shown in 
Fig.l). The value of diagonal exchange is considered to 
be Jd = 32.7meV which has been obtained by LDA-I-U 
calculations in Ref.( [Q). For K < Kc the first excited 
state is a spin triplet, for K > Kc however the first excited 
state is not a spin triplet, rather it has spin zero like the 
ground state and therefore it should be interpreted as a 
charge excitation. This leads to a very different scaling 
of the energy gap for K > Kc (Figj^b) while the case 
K < Kc is very similar to the model without diagonal 
exchange (Fig.^a). The scaling behaviour in Fig.^b for 
K<0.35 is roughly linear which behaves like the former 
spin excitations while for K>0.4 the scaling behaviour is 
completely different. 

J.„ (K) 



FIG. 8. Effective superexchange coupling versus K. It 
shows a steep decrease close to the transition point Kc=0.39. 
Moreover the effective exchange goes to zero in the completly 
zig-zag orderd phase (K^ cxo). 
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We have also considered the effect of dimerization (5) 
of the exchange coupUng on the energy gap of a single- 
ladder. This is of interest because it seems now clear 
that the charge disordered ladders in the low tempera- 
ture structure of NaV205 discussed in the introduction 
have a large superexchange dimerisation due to the shift 
of oxygen positions |4|. Therefore we consider the in- 
fluence of the dimerisation for K<Kc in the disordered 
regime. This means that the exchange constant (J) of 
the first term in Eq.(|^) is replaced by J(l ± 7) in an 
alternating way. The dimerization parameter varies as 
7=0.01, 0.05, 0.1, 0.2, . . ., 0.6. Obviously an energy gap 
opens by turning on 7 which does not scale to zero with 
increasing N (N — > cx)) as it did for the undimerized case 
7=0. In Fig^b we have plotted the log-log plot of en- 
ergy gap versus 7. For small value of 7 the energy gap 
scales as A ~ J'' where rj — 2.25 ± 0.25. 



(a) 





0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 
1/N 



(b) 




FIG. 9. Spin gap energy versus 1/N in the dimerized 
single-ladder for K=0.2<Kc and J=30 meV. (a) Data 
are plotted for different value of dimerization parameter 
(7 = 0.01, . . . ,0.6). All curves show up a finite gap in the 
N ^ 00 (for 7 > 0.05 it is more apparent), (b) The extrap- 
olated energy gap is plotted versus r in a log-log scale. For 
small value of 7 (7 < 0.1) we observe the scaling behaviour, 
A = 7". 



V. SUMMARY 



For the 2D Trellis lattice, it depends on the relative size 
of intra- and inter-ladder effective Coulomb interactions 
K,K' and the size of K with respect to the intra-rung 
hopping ta whether in-line or zig zag structures are real- 
ized. This agrees with LDA-I-U total energy calculations 
which show that these two structures are very close in 
energy with the zig-zag structure slightly lower Q. In 
addition however we have shown that in a wide region 
in the K, K' plane a disordered phase should exist which 
does not show long range correlations of charge due to the 
inherent frustration of the Trellis lattice structure. Ac- 
cording to Liidecke et al Q in the real NaV205 structure 
only every second ladder is charge ordered and the inter- 
vening ladders stay in the MV state. This might be an 
indication that NaV2 05 is not far from the regime with 
a disordered state, i.e. close to a quantum critical point. 
We have also studied the effect of spin coupling on the 
charge order transition. We have shown that Kc has a 
monotonic dependence on the exchange constant(J), i.e. 
for J=30 meV it leads to 15 percent increase of the ID 
ITF resuh. 

We have addressed the issue of excitation type in our 
charge-spin model. As long as the spin exchange along 
the ladder diagonal is absent the excitations are singlet 
to triplet. By adding a diagonal spin exchange on ladders 
we will observe charge excitations in the ordered phase 
while it is still spin type in the disordered phase. 

From the finite size scaling analysis of the singlet 
triplet excitation gap we obtain the dependence of the 
effective inter-rung exchange on the ladder as function of 
the degree of charge order, i.e. as function of K. A steep 
drop of Je/ / is observed around the critical value Kc for 
charge order. Eventually it will be reduced to zero in the 
fully CO regime in qualitative agreement with perturba- 
tion calculations [|4 
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The main results of our exact diagonalisation investi- 
gation for the present NaV205 model are the following: 
For the single ladder the ratio of K/ta controls the onset 
of charge order which is observed in the evolution of the 
pseudospin correlation function (Fig.^). This is similar 
as in the exact ID Ising model in transverse field results 
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